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Abstract. We introduce a small slit into a planar domain and study the re- 
sulting effect upon the eigenvalues of the Laplacian. In particular, wo show 
that as the length of the slit tends to zero, each real-analytic eigenvalue branch 
tends to an eigenvalue of the original domain. By combining this with our ear- 
lier work |HlrJdg07] , we obtain the following application: The generic multiply 
connected polygon has simple spectrum. 
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In this paper we study the following singular perturbation problem. Let fl C 
be a bounded open set having Lipschitz boundary. Remove from Q the horizontal 
slit, Et = {{x,y)\ y = yo, |a; — a^ol < t}, centered at the point {xo,yo) G ^- Standard 
analytic perturbation theory shows that the eigenvalues of the Laplacian on the slit 
domain \ S( depend analytically on t for t > 0. (See Theorem 14.21 ) The main 
result of our study is the following. 

Theorem 1.1. Let t i-^ Ef be a real- analytic eigenvalue branch of the Laplacian 
acting on L^(fl \ St) with either Dirichlet or Neumann boundary conditions. Then 
Et converges to an eigenvalue of the Laplacian acting on L'^{Vl) as t tends to 0. 

We remark that the convergence of real-analytic eigenbranches is much subtler 
than the convergence of ordered eigenvalues. For example, as t tends to zero, the 
A:*'' eigenvalue of the Laplacian on the rectangle [0,i] x [0, 1/t] tends to zero, but 
infinitely many analytic eigenvalue branches limit to infinity. Note that our proof 
that each eigenbranch in Theorem 11.11 has a finite limit is more involved than the 
proof that ordered eigenvalues limit to an eigenvalue of Q,. (Compare ^l3]and [HI) 

This work began as part of a study of the Laplace spectrum of degenerating 
translation surfaces. Indeed, the collapse of a slit is one of the typical degenerations 
of a translation surface |EMZ03| . Theorem II . II deals with the simplest case of such 
a degeneration. The slit degeneration of a translation surface is analogous to the 
well-studied 'neck pinching' degeneration of hyperbolic surfaces and our overall 
strategy for proving Theorem 11.11 mirrors the strategy employed in |Wlp92| and 
[Jdg02l . 

The strategy is to control the negative variation of the logarithmic derivative, 
dt\og{Et), for small t. A well-known formula — Proposition 14.61 — relates dtEt to 
certain quadratic forms defined as integrals over the domain. Not unexpectedly, 
one need only consider the contribution to these integrals that comes from a small 
neighborhood of the slit (CoroUarv 14. lip . A judicious choice of coordinates in a 
neighborhood of the slit allows one to reduce the analysis of the quadratic forms to 
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a family of 1-dimensional problems indexed by i = 0, 1, 2, 3 . . .. A simple convexity 
estimate (Lemma I5.3P provides for control of the 1-dimensional contributions for 
i > {2ty ■ Et where 2t is the width of the slit. We first show that ■ Et is bounded 
and then use convexity for i large and compactness for i small to find that t^^ ■ Et 
tends to zero for some fc < 1. Thus convexity applies to all z > 0, and a special 
estimate can be made for the case i — Q. 

Elliptical coordinates in a neighborhood of the slit are particularly well-suited 
for our purpose. Define (r, 9) implicitly by 



X — xo = Vr^+t^ • cos(0) 
y - yo = r ■ sm{9) 

The slit then corresponds exactly to the locus r ~ and the level sets r = const 
correspond to ellipses that surround the slit. If, in turn, one sets r = tsmh{z) then 
the equation Aip = E ■ ip \s separable in the variables (z, 9). The solutions to the 
resulting ordinary differential equations are called Mathieu functions. Appendix Rl 
contains the basic facts about these functions that we use here. We note that the 
work of Y. Colin de Verdiere [CdV87] called our attention to the usefulnes of these 
coordinates when considering slits, a fact that mathematical physicists have long 
been aware of [MrsRbn38] |MrsFsh| . 



We use Theorem 11.11 to extend the generic simplicity results of HlrJdgOJ ■ For 
example, we consider multiply connected polygons with n exterior verticesQ In §7, 
we prove the following: 

Theorem 1.2. //n > 4, then almost every multiply connected polygon has simple 
Laplace spectrum. 



In Theorem 1 1.21 we assume that the boundary conditions arc cither Neumann on 
every edge or Dirichlct on every edge. We also consider 'mixed' boundary conditions 
in the sense that the condition on each boundary segment is either Dirichlct or 
Neumann. 

Theorem 1.3. For each combinatorial choice of mixed Dirichlct- Neumann bound- 
ary conditions, almost every simply connected polygon with n > 4 vertices has 
simple Laplace spectrum. 

This latter theorem may be extended to multiply connected polygons with some 
mixed boundary conditions. It should be noted, however, that our method doesn't 
give generic simplicity for every type of mixed boundary conditions in the case of 
multiply connected polygons. For instance, we cannot handle a slit with Dirichlct 
conditions on one side and Neumann conditions on the other. 

We now outline the contents of this paper. In §2 we precisely describe the eigen- 
value problem on a slit domain and give a fuller description of elliptical coordinates. 
In §3 we prove that sequences of eigenfunctions with uniformly bounded eigenval- 
ues have convergent subsequences as the slit parameter t tends to zero. In §4 we 
apply Kato's theory of analytic perturbations to show that the eigenvalues and 
eigenfuntions depend analytically on t. We also derive variational formula for these 
eigenvalues and prove that uniform estimates on the logarithmic derivative can be 



^One can regard such a polygon as the result of removing finitely many simply connected 
polygons from a simply connected n-gon. 
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'localized' to a neighborhood of the slit. In §5 we prove a convexity result for 'ra- 
dial' Mathieu functions. In §6, we prove Theorem 11.11 according to the strategy 
describe above. In §7, we precisely state and then prove results concerning the 
spectral simplicity of polygons including Theorems 11.21 and 11.31 



2. The eigenvalue problem on a slit domain 

Let f2 C be a compact domain with Lipschitz boundary, and let S be a disjoint 
union of line segments contained in the interior of O. The complement O \ E is not 
a domain with Lipschitz boundary, and hence one must take care in defining the 
eigenvalue problem. The purpose of this section is to describe a resolution that will 
be useful to us. We refer the reader to [Grij for a complementary discussion of the 
analysis of singular domains. 

To define the eigenvalue problem, we will regard a slit domain as the compact 
manifold with Lipschitz boundary obtained by completing f2 \ E with respect to the 
length metric. To be more precise, we define the distance d{x, y) between x and y 
to be the infimum of lengths of rectifiable paths 7 : [0, 1] — > J7 \ E that join x to y. 
Note that d{x, y) equals the Euclidean distance between x and y iff a; and y belong 
to a convex subset of \ E. 

Definition 2.1. The slit domain fis is the metric completion of ri\ E with respect 
to d. We will refer to fJs \ (fJ \ E) as the slit. 

Although rij] is not isometric to a subdomain of , it is naturally a compact 
Ricmannian manifold with Lipschitz boundary. To see this, one can use elliptical 
coordinates to define a chart in a neighborhood of the slit. 

2.1. Elliptical coordinates in a neighborhood of the slit. Let = IR/(27rZ), 
and for each < > 0, define 6, x ^ 



4)t{r,9) = (^Vr2 +i2 .cos(6l), rsin(6l)) 



This injective map sends each circle {r} x onto an ellipse with foci at (±i,0). 
We have (/)i(M+ x S^) = \ where Ei = [-t,t] x {0}. By regarding \ Ej as 
a subset of Rl^, the map 0t extends to a smooth diffeomorphism from [0,oo) x 
onto R|^ . Abusing notation slightly, we will use 0^ to denote this extension. 

In these elliptical coordinates (r, 9) about the slit Et , the gradient operator takes 
the form 



„ , +t^ dw 1 dw 

^ ■ t^sin^{e) dr r^+t^sm^{9) 



and Lebesgue measure is expressed as 

(2) dm, = !^±i!£!!!!W drd9. 

The eigenvalue problem is unchanged if we rotate and/or translate the domain 
and the slit E C 51 simultaneously. Thus, in the sequel, we will often make the 
following assumption. 

Assumption 2.2. We assume that E = [— x {0}. 

Now choose ro > so that the ball B{0, 2ro) is contained in 51. Then for t < ro, 
the restriction of to [0, 2ro) x provides a chart in a neighborhood of the slit. 



4 



LUC HILLAIRET AND CHRIS JUDGE 



2.2. The eigenvalue problem. Let dm denote Lebesgue measure. For each 
smooth function on the manifold fis define 

N{u) = / \uf dm 
Jn 

and 




Given a measurable subset, D, of the boundary of Sis, define H\,{Q,y,) to be the 
completion of 

{u e C°°(f^E) I u{m) = 0, Vm e D} 
with respect to the norm 

u ^ q^{u) + N^{u). 

The form q extends to a closed quadratic form on H}-,[Q,y) and the form N extends 
to a closed quadratic form on L'^{U^,dm). In the sequel we will let n(-,-) (resp. 
q(-, •)) denote the polarization of N (resp. q). 

A function ip e Hjy{fls) is an eigenfunction with eigenvalue E if and only if 

q(?/;,u) = E ■ n{tp,v) 

for all V G iJ^(r2s). Integration by parts and standard elliptic estimates give that 
"0 is smooth in f2s with A-tp = E ■ where A is the Laplacian acting on smooth 
functions on M^. Moreover, we have 

(1) ip{m) ~ for all m G Int(_D) (Dirichlet conditions), and 

(2) v{4>){rn) = for all m S Int(9r2s \ D) (Neumann conditions). 

Here v denotes the outward normal derivative along the boundary, 9f2s, of fis. 
Note that in the sequel D will essentially be a union of segments. 

3. Convergence of eigenfunctions with bounded eigenvalues 

Let r2 be a domain that contains the origin {0}, and let t„ > be a sequence 
with lim„^oc tn = 0. For sufficiently large n, the segment Yit^ — [~tn, i„] x {0} lies 
in the interior of il. For each n e N, let Dn be a measurable subset of the boundary 
of the slit domain ^Is^^ such that D = £)„ n d^l does not depend on n. Let ipn be 
a normalized eigenfunction of q on Hj^^i^'St^ ) with eigenvalue En- In this section 
we prove the following: 

Theorem 3.1. If E„ is a bounded sequence, then a subsequence ipn{k) converges 
in L'^{fl) to an eigenfunction ip of q on Hj^(n) with eigenvalue E ~ lim^^oo ^'rt(fc)- 
Moreover, for every neighborhood U of the slit and every j e N, the sequence ipn{k) 
converges to -0 in (fl \ U). 

Proof. Let J7 be a neighborhood of the slit. Each ipn is an eigenfunction of A with 
eigenvalue En, and hence for each n such that C C/ we have 

/ I^Z-n • A-'V'nl dm = E^n / \ipn\'^ dm. 
Jn\u Jn\u 

Thus, using Garding's inequality and the fact that 0„ is normalized, we find that 

ll^nll, < c.(£;^ + i) 
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where || • \\j is the norm associated to the Sobolev space H^{Q, \ U). Since by 
assumption, En is a bounded sequence, we have that for each j, the sequence 
n ^ llV'nIlj is bounded. 

Using a diagonahzation argument, we find a function -0 and a subsequence of 
■0n — still denoted -i/jn — such that for every j and every neighborhood U of the slit, 
ipn converges to i}} in the (Cl\U) norm. Thus, by the Sobolev embedding theorem, 
tpn converges to tp in C''{il\U) for every k £ N. 

By assumption, the ipn have L^-norm equal to 1. Thus, given e > 0, Lemma 13.21 
below applies to give A^i > and r* so that if n > A^i, then 



where B is some fixed ball centered at and included in the neighborhood Ut. 
for n > Ni. Hence, by Fatou's Lemma, we also have 



Jb " 

The functions ipn converge uniformly to ijj on the complement of B, and hence there 
exists N2 so that if n > iV2 then 



It follows that V'ri converges to ip in L^(f2). In particular, ?/; is i^-normalized and 
hence is nontrivial. 

Finally, we show that ip is an eigenfunction. To do this we adapt an argument 
from |CdV82| . Define the distribution T by T{(j>) = J^ip ■ (p dm. Then A'T e 
H^^iyi) where A' denotes the distributional Laplacian. If the support of (j) does 
not contain the origin, then A'T{(j)) = E ■ T(0), and hence the singular support of 
A'T-ET is contained in {0}. Hence there exists S € H~^{n) with supp(S') C {0} 
such that 



It follows that there exists C such that S — C ■ 5 where = 0(0). 

Let G be the distribution defined by G{(j)) = J^L ■ cj) dm where L{x,y) = 
ln{^yx^ + y^)x{x,y) where x is some smooth cut-off function near (say for in- 
stance in the ball B). A direct computation shows that A'G — (5 is in L^(il) and 
hence it follows from above that A'(T — C-G) is also in L'^{fl) and hence, by elliptic 
regularity, T - C ■ G € H^{Vl). 

It suffices to show that (7 = 0. On the one hand, by Fatou's Lemma, we have 






Q.\U 



A'T ^ E-T + S. 




= E. 



and on the other hand 



/ |VLp dm > 2tt f 



{dr ln(r))^ r dr 



= 00 




where -8(0, e) C is a ball centered at the origin. Therefore, T — G ■ G ^ H^{il) 
unless C = 0. □ 
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Lemma 3.2. Given e > and Eq > 0, there exists r* > such that if u is an 
eigenfunction of q on Hj~)(fl^^) with eigenvalue E < Eq, then for allt < r* we have 

|mP dm < e / |mP dm 

where Ut^r' is the elliptical neighborhood of the slit T.t of radius r* . 
Proof. For positive r and p, we have 

\u{r + p,e)-u{r,e)\ < f \dru{r + s,9)\ ds 



< Vp- (^J^^ \dru{r + s,9)\^ ds 
by the Cauchy-Sehwarz inequaUty. From this we find that 

(3) l-\u{r,9)\' < \u{r + p,9)\' + p j'^ \dru{r + s,9f ds. 

Given ri > 0, define Ug — {r \ s < r < ri + s} x . Using ^ and @, we find 
that if s > 

\dru{r + s,9)\^ dm < [ \dru{r + s,9)\^ {r^ + t^)^ d9dr 

< [ \dru{r + s, 9)\^ ((r + sf + t^)^ d9dr 

Jua 

|9,M(r,0)|2 (r2+t2)^ d9dr. 

< I \Wu{r,9)\'^ dm. 

JUs 

Thus, since u is an eigenfunction with eigenvalue E, 

(4) P ( \dru{r + s,9)\^ dmds < p^E [ \u\^ dm. 

Jo Juo Jn 

Fix some p so that p^E < e/4, and choose r* < p. Uniformly for r,t < r* , we 
have then 



+ sin^ ' 



'-T^,< (r^ + t^)^ < r* 
(r2 +i2)^ - ^ ^ - 

and 

(r + p)^ + sin^ 9 (r + p)^ p 

(r + /3)2+t2)i - (2(r + p)2)i - 
from which it follows that 

r2 + f2 sin^ 9 ^ r* \/2 (r + p)2 + i2 sin2 Q 



(r2 + i2)^ - p i^i^r + p)^ + t^)i 



We now fix r* < so that we have 



r2 + <2 sin2 9 ^ e {r + pY + t^ sin^ ( 



uniformly for t <r* . 
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Thus, using ^ and (O, we have 



Uo 



, ,u(r,d)\ ■ — dOdr 



< ^ / |u(r,0)|2 dm. 



By combining this with Q and ([4]) wc find that 

The claim foUows since p^E < e/A. □ 



4. AnALYTICITY, VARIATIONAL FORMULAE, AND LOCALIZATION 

In this section we use standard perturbation theory to show that the eigenvalues 
and eigenfunctions depend analytically on t for t > Q. We then derive basic formulae 
relating the derivative of an eigenbranch to the derivatives of the quadratic forms q 
and N . We show that for small t the dominant terms in the derivative fomulae can 
be localized to a neighborhood of the slit. Finally, we evaluate these local formulae 
on an elliptical neighborhood of the slit. 

4.1. Pulling back to elliptical coordinates. Standard analytic perturbation 
theory [Katoj applies to a family of quadratic forms on a fixed Hilbert space. For 
this reason we will modify the family quadratic forms q and N but in such a way 
that we obtain an equivalent eigenvalue problem. 

Of the various possible approaches, wc choose to modify the map (pt so that the 
inverse image of is constant in t. The fixed Hilbert space will then be on 
this inverse image with respect to the pull-back of Lebesgue measure. 

Let To be as in fJ2] Namely, -6(0, 2ro) C 51. Let tr : R ^ M be a smooth, positive, 
decreasing function such that a{t) = 1 if |t| < 1 and a{t) = if \t\ > 2. Define 
0t : R X -> R2 by 

4it{r,9) = (^y^r^ + ■ (j{r/ra) ■ cos{9), rsm{9)j . 

By construction, (pt defines a smooth diffcomorphism from [0, oo) x R onto R|;^. 

If i < 7'o, then the set M = does not depend on t. We pull-back the 

Dirichlet energy functional and the L^-norm on to functionals on M. 

For computational convenience, we express the change of variables in the lan- 
guage of Riemannian geometry. For each t, define the Riemannian metric 

gt = ^;{dx^ + dy'') 

on [0, oo) X . In particular, is a Riemannian isometry from (fit, dx^ + dy'^) onto 
iM,gt). 

Let dmt denote the Riemannian measure: dmt = {,4't)Z^ {dxdy) . The space, 
L'^{M), of functions that are square integrable with respect to dmt does not depend 



8 



LUC HILLAIRET AND CHRIS JUDGE 



on t. On the other hand, the associated norm does depend on t. Let Nt denote the 
quadratic fornJl 



(5) N,{u) 

J M 

In the sequel we will let nt(-, •) denote the polarization of Nt- 

Let Vt denote the Riemannian gradient for (M, gt). Namely, we have gt{WtU, X) = 
du{X) for every smooth function u : M -^M. and vector field X on M . For smooth 
functions u : M ^M. define 

(6) qt{u) = / gt{VtU,Vtu) dmt. 

J M 

Given a measurable subset D in the boundary of Af, define H]j{M) to be the 
completion with respect to the norm 

of the set of smooth functions on M that vanish on D. Since the equivalence class 
of this norm is independent of t > 0, the completion does not depend on t. The 
form qt extends to a closed quadratic form on H]-,{M). In the sequel we will let 
qi(-, •) denote the polarization of qt- 

Proposition 4.1. ijj is an eigenfunction of the quadratic form q on i7^^j^j(f2st) 

with respect to N on L^(rist) with eigenvalue E if and only if u ~ 4>t{i') '•s 
eigenfunction of qt on H}j{M) with respect to Nt on L^{M) with eigenvalue E. 

Proof. This follows from a straightforward accounting of the change of variables 
given by 0t. □ 

4.2. Analyticity of eigenvalues and eigenfunctions. 

Proposition 4.2. The eigenvalues and eigenfunctions of qt with respect to Nt vary 
real- analytically. To he precise: V/c e N, 3 real-analytic paths i/jk ■ (0,to] Hjy{AI) 
and Ek '. (0, to] ^ ^ such that 

(a) ([t{Mt).v)=Ek{t)- nt ( V-fe (t) , for all v e Hj, (A/) . 

(b) NtiMt)) = 1. and 

(c) the span of {tpkit) \ k £ N} is dense in L'^{M). 

Proof. Apply standard analytic perturbation theory. In particular, both qt and 
Nt are holomorphic families of quadratic forms of type (a) as in §VII.4 |Katoj . 
The eigenvalue problem is of generalized form as in §VII.6 jKato| . (See especially 
Remark VII. 6. 2 |Kato| and the discussion on page 419.) The analogue of Remark 
VII.4.22 [Kato] gives the claim. □ 

Remark 4.3. We will use the expression eigenbranch to designate the mapping 
t (ipkit), Ek(t)) and will use the expression normalized eigenbranch if, in addition, 
Nt{ipk{t)) = 1. By extension, the term eigenbranch may also refer to either the 
eigenvalue or the eigenvector singly. Since we will be dealing with one eigenbranch 
at a time, the index k will be systematically dropped. 



"^This quadratic form is defined on functions on M, and henee is not the same as the quadratic 
form A'^ defined in [[2] On the other hand, these forms differ by pull-back by tpt, and the context 
will make clear which is being used. 
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Corollary 4.4. The Neumann (resp. Dirichlet) eigenvalues of the Laplacian on 
f2Et depend analytically on t G (0,io]- 

Proof. Combine Proposition [01 Proposition [521 a-nd the discussion in 321 O 
4.3. Variational formulae. 

Notation 4.5. A dot above a quantity will indicate differentiation with respect 
to t. For example, E indicates the first derivative of an eigenvalue branch Et. In 
what follows, we will often suppress the dependence of q, ipi ^^nd E on t from the 
notation. 

We begin with a well-known, general variational formula. 
Proposition 4.6. We have 

E-N{%P) = q{i') - E-NitP). 

Proof. Substitution of i/' for v in part (a) of Proposition l4. 21 gives 

(7) q(i^,^) = i?-n(V',V'). 
By differentiating part (a) of Proposition [521 we obtain 

q(i/',w) + q(V',w) = £'-n(V', w) + E ■ {n{i\v) + n(?/', w)^ 
for all V. By substituting ip for v and using ([7|), we find that 
qiip.ii) = E-N{iP) + E-N{iP). 

□ 

In the case that q is the Dirichlet energy associated to a family of Riemannian 
metrics, the quantities in Proposition 14.61 can be expressed in terms of the first 
variation of the metric and the associated Riemannian measure. 

Proposition 4.7. Let t —f gt be a real-analytic family of Riemannian metrics on 
M , and suppose that q is defined by (0). Then we have 

q{u) = ^ g (Vu, Vii) dm + g (Vm, Vm) dm 
JM Jm ' 

and 

N{u) = [ |up dm. 
Jm 

Proof. By differentiating ^ we have 

q{u) = / g (Vu, Vu) dm + 2 / glVu, Vuj dm + / 51 (Vu, Vm) dm. 
Jm Jm ^ ' Jm 

By definition g(Vu,X) = ^{f) for all vector fields X, and hence we have 

(8) .9(V«,X) + g(Vu,X) = 0. 
In particular, if X = Vu, then we have 

/ g{Vu,Vu) dm + / g{Vu,Vu) dm = 0. 
Jm Jm 

Substitution into the formula for q gives the first formula. The second formula 

follows from ([5|). □ 
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The following lemma will be used to translate estimates on the logarithmic de- 
rivative of an eigenbranch into statements concerning the convergence of the eigen- 
branch. 

Lemma 4.8. // there exists to and a continuous positive function p : (0, to] 
such that 

dtEt > -p{t)-Et 
for all t S (O,to]i then the function 

rto 



Fit) = exp {- pis) ds^ ■ Et 



converges as t tends to zero. In particular, if p is integrable, then Et converges as 
t tends to zero. 



Proof. We have 

pis) ds ] > 0, 



F'it) - idtEt + pit) ■ E) ■ exp 

and hence F is increasing. Since F is nonnegative. the claim follows. □ 
Using this lemma and Proposition 14.61 we have the following corollary. 



Corollary 4.9. Suppose that there exists a constant C such that for allu £ H^iM) 
we have 

(9) qiu) > -C-qiu) 
and 

(10) Niu) < C-Niu). 
Then Et converges as t tends to 0. 



Proof. From Proposition |4]6] we have E > — 2C • E. Thus we can apply Lemma 14781 
with p = 2C. □ 

4.4. A localization principle. Unfortunately, inequalities © and pI7|) do not 
hold true for the singular perturbation that we consider here. We now show, how- 
ever, that these inequalities do hold for all u with support outside of a neighborhood 
of the slit (Proposition HUO]). 

To state this result in a convenient form, we introduce the following notation. 
Let [/ be a measurable setH For w G H^, we denote by quiw) the 'restriction' of q 
to U. That is, 

quiw) = / iVwpdm. 
Ju 

We will use analogous notation for the quadratic forms q, N , and N . 



Proposition 4.10. Let U G be a neighborhood of the slit. There exists a constant 
Cjj such that for any w € 

qM\uiw) > -Cu ■ qiw) 

and 

NM\uiw) < Cu-Niw). 



•^In the sequel, U will often be an elliptical neighborhood of the slit. 
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Proof. Let SM denote the unit tangent bundle to M (with respect to gt,, ) , and for 
each smah t, define Ft : SM ^ R by 



Ft{X) 



g(x,x) 

The restriction of gt to il\U is real-analytic for t E [0,to], and hence 

Ci = sup{Ft{X) I t G [0,io] and X e S{n\U)} 

is finite. By homogeneity of g and g in each tangent space, we have g{X,X) < 
Ci ■ g{X, X) for all X € T{n \ U). Therefore, 

/ giyw^Vw) dm < Ci g{\7w,\7w) dm 

JM\U Jm\u' 

Let Gt : [0, r^] x 71/ ^ M be defined by 

Since dmt restricted to O \ [/ depends real-analytically on t e [0, to], 

Ca = sup{|/(t,X)| I t e [0,ro] and X en\U} 
is finite, ft follows that 



and 



ui^ dm < 6*2 / dm, 

/ g (Vw, Vui) dm > — C2 / g (Vw, Vw) dm. 
Jm\u' Jm\u 



lM\U JM\U 

We use Proposition l4.7l and set Cjj — Ci + C2- The result then follows. □ 

Corollary 4.11. Let U be a neighborhood of the slit. Suppose that there exist 
functions a,(3,j : [0,oo) — > [0, 00) such that the following holds for an eigenbranch 

quii^) > -a{t)-quii^) - p{t) ■ E ■ NuO^), 

and 

NuW < i{t)-NuW. 
Then there exists C" > such that 

E > -{C - a{t) - Pit) -i{t))E. 

Proof. We have 

= 9a/\c/(^) + 9c/ (-0) 

> -{C + a)-q{i^) - p-E-N{ij) 

> -{C + a + P) ■ E ■ N{i') 

where the constant C comes from Proposition 14. 101 Similarly, we find that 

The estimate then follows from Proposition 14.61 with C = 2C. □ 
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Remark 4.12. Proposition l4.10l and Corollary [47TT] represent a localization principle 
in the sense that an estimate of the logarithmic derivative of Et has been reduced 
to the study of the functional q and in a small neighborhood of the slit. This 
principle applies more generally to singular perturbation problems in which the 
'singular support' of the perturbation is small. 

4.5. Evaluation on the ellipse. Wc now let U be the elliptical neighborhood 
of the slit of radius ro and evaluate qu and Nij. By using the expression for the 
gradient H]) and the Lebesgue measure ([2]) in elliptical coordinates, we find that 

q^(v) ^ f \drv\'ir'+t')'^drd0 + f \dovf 

Ju Ju [r^ +1-^)2 

and 

f , ,Jr2+t2sin2(0)) 

Ju {r^+t^)2 

By Proposition 14. 71 or direct computation, wc find that 

drd0 /" , „ ,1 drdO 

quiv) = t / \drv\' — -X - t / \d0v\ 



and 



2((2sin2(6l)-l)r2+t2 sin2(e')) 



Nu{v) = t / \v\'^ ^-^^^^ ^ drde. 



It will prove convenient to make a change of variables. 

Notation 4.13. In the following, wc let 
m r ^ t ■ sinha;, 

• Yt ^ sinh^^(ro/t), and 

• Ut = [0,Yt] X S\ 

With this change of coordinates, the formulae above become 

(11) qu{v) = [ dxde + [ \dgv\^ dxde, 

Jut Jut 



(12) Nuiy) = i 



[ \v\^ {sinh^ x + sin^ {9)) dxde, 
Jut 



(13) qu{v) = i / ^ dxde 1 / ^ dxde, 

t Jut cosh x t Jut cosh x 

and 

(14) Nu{v) = t I \v\^{&\\?{e) -cos2(6l)tanh2(.T)) dxdO. 

Jut 

From HH) and wc find that 

(15) quiv) > - 7 / dxde, 

t Jut cosh X 

and 

(16) Nu{v) < t [ \v\^sm^ (9) dxde < t [ \v\^ dxde 

Jut JUt 
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5. MATHIEU functions and INTEGRAL ESTIMATES 

In this section wc provide estimates that wiU be used to control ratios such as 
N /N and q/q in the next section. 

5.1. Separation of variables. E. Mathieu [Mth68| observed that one can apply 
the method of separation of variables to the eigenvalue problem, A^/; = • -0, on an 
ellipse in . He made a detailed study of the solutions to the resulting ordinary 
differential equations, solutions that now bear his name. To perform the separation 
of variables, we use the same change of variables as in the preceding section, i.e. 
we set 

r = t ■ sinh(a;). 

Let /i £ M. The operator 

acting on L^(M/(27rZ), rffl) has discrete spectrum 

< bf){h) < bi{h) < b2{h) < ■■■ . 

For h ^ 0, the spectrum is simple [MrsFsh] . Let Vi^h denote the i^-normalized 
eigenfunction associated to bi{h). We will call Vi^h the i*^ angular Mathieu function. 
In the following we will often suppress the dependence of Vi^h on h from the notation. 
For the convenience of the reader, we prove the following in Appendix [K\ 

Proposition 5.1. Suppose that Atjj = E-tp on the ellipse of radius rp = t sinh(a;o), 
and suppose that tp satisfies a Dirichlet (resp. Neumann) boundary condition on 
the slit: For Q < 6 < 2-k we have 

^-(0, 6*) = {resp. d^ipiQ, 9) = 0). 

Then for x > 0, 

(17) ^|;ix,e) = ^ii,(x) •z;,(0) 

i 

where 

h ^ t - /e, 
and where ui : [0,2:o) ^ M is a solution to 

(18) - u"{x) + {b,{h) - h^ coiih^{x)) - nix) = 
with w(0) = (resp. it'(0) Oj. 

A solution u to (|18p will be called a radial Mathieu function. 

In the sequel, our methods will rely upon estimates of the distribution of the 
and mass of an eigenfunction on the ellipse. The following lemma will allow us 
to reduce such estimates to estimates of radial Mathieu functions. 

Lemma 5.2. We have 

(19) / de 

(20) / \deijfde + IVP cos\0) dO ^ bdh) ■ \u,{x)\'' . 



14 



LUC HILLAIRET AND CHRIS JUDGE 



Proof. By standard Sturm-Liouvillc theory, the angular Mathieu functions {vi} 
form a complete and orthogonal set in L?{S^,d9). (See Appendix E}. Thus, HH) 
follows from P7|) . 

From (I17p . we have dg^ — "^Ui-dgVi. Since Vi is an eigcnf unction with eigenvalue 
bi = bi{h) we have 

-d^Vi + COS^(0) • = fei • W,;. 

By multiplying by Vj and integrating by parts, we have 

dgVi ■ dgVj d9 + h'^ / Vi ■ vj coa'^{6) dO = bi j Vi ■ vj dO. 

Equation (jSO]) then follows from the fact that {vi} is complete and orthogonal. □ 

5.2. Radial convexity estimates. Our estimates of the distribution of the 
mass of a radial Mathieu function u depend upon the following observation. We 
have 

(u^)" = 2-u-u" + 2-iu'f > 2-u-u" , 
and hence if hi — cosh^ (x) > ^ , then by we have 

> u\ 

We will let w : [0, X] M+ denote a smooth function such that 

w" {x) > 'w{x), 
w{x) > 



for aU X e [0,X] and 



y'(0) > 0. 



In particular, if X satisfies 

(21) 6, - •cosh^(X) > i, 

then the square, uf, of a radial Mathieu function satifying either Dirichlet or Neu- 
mann conditions is an example of such a function w. 

The following expression of convexity is the basis for our estimates. 

Lemma 5.3. For all x,y > such that x + y < X , we have 

w(x + y) > w{x) - cosh{y). 

Proof. The claim holds if w{x) = 0. So we may assume that w{x) > 0. Let z{y) = 
w{x + y)/w{x). Note that z{0) = 1 and z'{y) > 0. Since w"{x + y) > w{x + y), we 

have -^(^(j/) — cosh(y)) > and since w'(x) > 0, we have 4- (z(y)~coah(y)) > 
ay ay 

0. It follows that dy{z{y) — cosh(j/)) > for all y >Q. Since z(0) — cosh(O) = 0, we 

have z{y) — cosh(2/) > for all y, and the result follows. □ 

Proposition 5.4. Let p : [0,X] —t M+ he a decreasing integrable function. Then 
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Proof. Applying Lemma 15.31 with y = X/2 gives 

cosh(X/2) / w{x) dx = w{x + X/2) dx < / w{x) dy 

Jo Jo J x/2 

and hence 

fX/2 ^ j-X 



Since p is decreasing 

X pX/2 pX 

p{x) ■ w{x) dx < p(0) / w{x) dx + p{X/2) / w{x) dx. 
Jo J x/2 

Combining this with ((22|) gives the claim. □ 



Proposition 5.5. Let p : [0,X] —t M+ he an increasing integrable function. Th- 



en 



X 2 rX 







w{x) dx < / w{x) ■ p{x) dx. 

p(X/2) 7o 



Proof. From (f22|) we have 



w{x) dx < ( rTTTTTTT + l) / w{x) dx. 

Vcosh(X/2) J J x/2 



Since p is increasing, we have 



X ^ pX 



w{x) dx < / w{x)-p{x)dx. 

x/2 P{XI2) J x/2 

By combining these inequalities and using the fact that cosh(X/2) > 1, we obtain 
the claim. □ 

6. Limits for analytic eigenbranches 

In this section we prove that each real-analytic eigenvalue branch Et converges. 
The proof consists of three main steps. First we prove that t^ ■ Et converges as 
t tends to zero. We use this to then prove that t^^Et converges for some k < 1. 
Finally, we use this to prove that Et converges. At each new stage, the previous 
estimate is used to control h{t) along the eigenbranch. 

6.1. Convergence relative to t^. 

Proposition 6.1. Let Et he any eigenhranch. Then 

lim ■ Et 

t^o+ 

exists and is finite. 

Proof. Since cosh^(a;) > 1, by comparing (lllll and (fT5|l we find that 

t-quiip) > -q{-ip) 
and by comparing (fT2|) and (fT6|) 
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Hence by Corollary 14. 1 II we have 

(23) E > - (^C + j^E 

for some constant C. Lemma l478l then allows to conclude since m this case F{t) = 
Ct^Et for some constant C . □ 



6.2. Convergence relative to t^^ . By Proposition [611] the parameter h ~ t - \fE 
in the radial Mathieu equation associated to -0 is uniformly bounded in t. This will 
allow us to prove the following. 



Theorem 6.2. Let Et be any eigenbranch. Then there exists k < 1 such that 



hm t^^ ■ Ef 



exists and is finite. 

Proof. It suffices to show that there exists k < 1 so that 

(24) t-NuW < n-N{i})). 

For then we could argue as in the proof of Proposition lOl where the '2' that appears 
in ([251) is replaced by k + 1. In particular, the desired k equals (k + l)/2. 

Lemma 6.3. There exists M > such that 

I dxdO < M I IV^p sinh^ xdxd0. 

Jut Jut 

Assuming the lemma, we finish the proof of the theorem. Let < k < 1 be such 
that M = k/(1 - k). Then 

(1-k) / ItPl^ sin^ 9 dxdd < (1 - k) / \tP\^ dxd9 < k [ li^l^ sinh^ x dxdO 
Jut Jut Jut 

and hence 

/ sin^ 9 dxdO < K I sinh^ X dxd9 + k [ sin^ 9 dxd9 . 
Jut Jut JUt 

Estimate ([24|l follows then by comparing ([T2|) and ^6]) . □ 

Proof of Lemma\KM By Lemma [5^ it suffices to prove that there exists M so that 
for all i 

(25) / \ui\'^ dx < M / sinh^xdx. 



Let p{x) = sinh^(x), and let X = sinh~^(l). Since p is increasing on [0,oo), the 
infimum of p{x) over [X, oo) equals p{X) = 1. In particular, if i < tq, then Yt > X, 
and we obtain 

(26) / \u,f p{x) dx > \u^f dx. 

Jx Jx 

By Proposition 16. li h ~ h{t) is uniformly bounded. Choose io G Z+ so that 
*o ^ 2h'^{t) + 1/2 for all sufficiently small t. Then since bi{t) > i^ — see Appendix 
[Al-and cosh2(X) = 2, we have that for all x E [0, X] and i > io 

b^{t) - h{tf ■ cosh^(x) > i. 
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Since '0 satisfies either Neumann or Dirichlct conditions along the slit, we have 
either u';(0) = or Ui{Q) = for all i. Note also that p is increasing. 

Thus, for i > ioj we may apply Proposition 15 . 51 with w = uf to find that 

The lemma is thus proved for i > io large enough. 

For i < io, we note that by Proposition 16.11 h is bounded, and hence bi is 
bounded. Thus, the claim for i < io follows from the fact that the solution of 
the ordinary differential equation (|18p with a fixed boundary condition depends 
continuously on parameters. □ 

In the next section we make crucial use of the following variant of Theorem 16.21 
Corollary 6.4. Let h{t) = t ■ s/E't. There exists eg > so that 

lim = 0. 

6.3. Convergence of Ef. 

Theorem 6.5. Let Et he a real- analytic eigenbranch. Then 

lim Et 

exists and is finite. 

Proof. It suffices to show that there exists (5 > so that for all sufficiently small t 

(27) quW > - t'-' ■ quW - 2t'-' ■ E ■ NuW, 
and 

(28) NuW < t^-^-Nui^P). 

For then by Corollarv l4.11l there would exist C > so that 

(29) E > - {C + A-t^-^)E. 



Since t^~^ is integrable near t = 0, the claim would then follow from Lemma 
By comparing pT|) to and to pB)) . we see that to prove (^7)) and 
it suffices to prove that there exists (5 > so that for sufficiently small t we have 

(30) 



and 
(31) 



/ dxde < t^ ( [ dxde + 2t'^ ■ E ■ [ ip'^ sinh^ x dxdo) 

Jut cosh X \JUt Jut J 

I \'il}t\'^ dxd9 < t^ I^Z-tp smh^x dxde. 
Jut Jut 

These estimates can be further reduced to estimates of Mathieu functions as in 
§5.21 In order to state these estimates, we first note that, by Corollary 16.41 there 
exists e > such that 

(32) ' ^'^^^ = ^• 

We fix this e in what follows. 
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Lemma 6.6. Let Ui be a solution to the radial Mathieu equation il8]} with h = hit) 
satisfying \32(l . There exists to so that if t < to, then for all i > 



(33) / \u,\^ dx < fl"^ I \ui\^smh^ X dx. 

and for z > 



(34) / dx < / \u,\'' dx 



cosh^ x 

Assuming this lemma, wc finish the proof of Tlieorem 16.51 Using (|19p we see 
that (|3ip foUows immediately from ((33l) with 5 = e/2. 
To verify ^ we use ^ to find that 

Jut cosh (x) ~^ Jo cosh (x) 

In Appendix [Xl we show that bo ~ i/i^ for small h. Hence since = t^E, we find 
that for small t 

50 < ^2 ^2(^) 

Jo cosh (a;) Jo 

rYt 

< E-t^ ■ t"/^ / ul{x) sinh^ a; dx, 

Jo 

< -t^ -E 1 IV'P sinh^ X dxde. 

Jut 

Here we have used (|33)) . Using ((20| and ([34]), we find that 

E / ^iS^ ^ / ^'(^) 

^ Jo cosh {x) ^ Jo 

< (^j \dgil;\'^dxde + h^ j \il;\'^dxde^ . 

Thus, since h"^ = t^E, and using to estimate the last integral of the previous 
inequality, equation ([50)) follows with S ^ e/2. □ 

6.4. The individual estimates. The remainder of this section is dedicated to 
proving Lemma 16.61 The behavior of the radial Mathieu function m for small h is 
ciualitativcly different if i = than if i 7^ 0. We prove the cases i ^ with Lemma 
16.71 and the case i — with Lemma [6791 

6.4.1. Estimates for i ^ 0. 

Lemma 6.7. There exists some to and some constants Ci,C2 such that for any 
t < to and for any i ^ 0, we have 

(35) / * dx < Ci-f r uj{x) dx, 

Jo cosh [x) Jo 



an- 



d 



(36) / u1{x) dx < C2-f u'^{x)smh^{x) dx. 

Jo Jo 
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Proof. Define Xt G [0, oo) implicitly by cosh(Xt) = t ^. Then we have 

d:c < t^' r ufix) dx. 
Jxt cosh (x) Jxt 

Since = o(t^). It remains to estimate the integral over [0,Xt]. 

By ((32)) . we have that h{t)'^ ■ cosh(Xi)^ tends to zero as t tends to zero. Since 
i ^ 0, we also have 6^ > 1. (See Appendix [X)) . Thus, for sufficiently small t, 

h{h{t)) ~h{tf -cosl^iXt) > i. 

Therefore, we can apply the estimates of ^5.'2\ with w = and X = Xt. By 
applying Proposition 15.41 with p = cosh~^, we find that 

dx < / uf{x) dx. 



cosh^(x) cosh(Xt/2) Jq 



Since cosh(Xt/2) = y^{cosh{Xt) + l)/2 = y^it-^ + l)/2, we have 

dx < — / uf{x) dx. 



cosh-'(a;) ~ I + 

Estimate ([55)1 then follows for Ci = 2\/2 and any i small enough so that < Ci-f^ 
The proof of the estimate ([36]) is similar. From the definition of Xt we have 



/ uf dx < — / uf ■ sinh(x)^ dx, 



and from CoroUarv 15.41 we find that 

Xt g^e rXt 

uf{x) dx < / sinh^(x) ■uf{x) dx, 

^ — Jo 

the estimate follows with any C2 > 8 and t small enough so that j^^pr and 

are less than C2. □ 

Since equation (|35p implies equation (j34p . and equation p6p implies equation 
(|33p . the estimates of Lemma 16.61 are proven for any index i > 0. It remains to 
prove ([55)1 for i = 0. 

6.4.2. The estimate for i — 0. Since 6o(^) is real-analytic in h^, we can write 

bo{h) = a{h) ■ 

for some analytic function a. It follows from Theorem 16.21 that a{h{t)) is bounded 
near t = 0. The Mathieu equation for i = can thus be rewritten as 

(37) u^{x) = h{tf ■ {a{h{t)) ~ cosh2(a;)) • ^0(2^)- 

We will treat as this equation as a perturbation of Uq — 0. To make this precise, 
we set io > and define 

M{X) = sup { \a{h{t)) - cosh^(a;)| | < a; < X and < t < io} • 
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Lemma 6.8. Given a solution u to \3'7\j define R by 

R{x) = u{x) ~ (u(0) + u'(0) • x) . 
Then for all t G (0,to), and for any positive weight p, the function R satisfies 



X 



X 



whe 



\R{x)\^ p{x) dx < h{tf ■ M{XY ■ ■ L{X) ■ / \u(x)Y dx 



L{X) = sup p{x). 

a:e[0,X] 



Proof. By applying the method of variation of constants (or, cquivalcntly, Duhamel's 
principle) wc find that 

R[x) = h[tf I [x~v) {a{h{t)) - cosh2(y)) u{y) dy. 
Jo 

Thus, for t <to and < x < X 

< h'-M{X)^(^J\x~y)-\u{y)\dy^ . 
By applying the Bunyakovsky-Cauchy-Schwarz inequality wc find that 



\R{x)\^ < h^-M^(X) 



{x - yfdy 



u{yf dy 



< h^-M^(X)-{\x\y3) 



u{yf dy 



for all |a;| < X. Thus, 

rX pX 

\R{x)\^ p{x) dx < h'^ ■ M^{X) ■ X^ ■ L{X) 



u{yf dy 



dx. 



The desired estimate follows. □ 

The following completes the proof of Lemma 16.61 

Lemma 6.9. Let e > Q be as in iSS^) . There exists some constant C such that, if 
we let u ~ uq be a solution to the radial Mathieu equation I118\} associated to bo{t) 
and h{t), and if u satisfies either u{0) = or it'(O) = 0, then for any sufficiently 
small t we have 



< C-t" 



Yt 



u'^{x) sinh^(a;) dx. 



Proof. As above, let > be defined by cosh(Xt) = t~'^. Let a — u'{Q) and 
h = u{0). Since < 2{ax + b)'^ + 2R{x)'^ , by applying Lemma with v(x) = 1, 

we find that 

Xt 

(38) 
where 

In other words. 



u^{x) dx < 2 {ax + by dx + K{t) / u\x) dx 



K{t) = 2M\Xt)-h{t)^-Xf. 



(39) 



Xt 



u'^{x) dx < 



1 - Kit) Jo 



Xt 



{ax + bY dx. 
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provided that K{t) is less than 1. As t tends to zero, Xt ^ e\ \nt\ tends to infinity, 
and hence - cosh^{Xt) - t'^". It follows that K{t) - (t"^ • /i • e| In 

Thus, by ([5^ -^(0 tends to zero as t tends to zero. Therefore, for t sufficiently 
small 

(40) / u'^ix) dx < 4 / {ax + b)'^dx. 



If w(0) 0, then ax + b = u'{0) ■ x, and if u'(0) = 0, then ax + b = u{0). A 
straightforward calculation gives a constant C such that for all sufficiently large X, 
we have 

f^X /~i \^ fX 



f C ■ X f 

i 1 dx < 2 / sinh^(a;) da:. 

Jo sinh {X) Jo 



and 



C ■ X 

/ x^ dx < ^ / sinh^(x) dec. 

Jo ~ smh^(X)Jo 



sinh^(X) Jo 

Thus, since linit^o Xt = oo, for small t we have 



Xt c V 

\2 ^ O • / ,^ , ,x2 2, 



(41) / {ax + bY dx < ^— — / (ax + bV sinh^x) dx. 

' smh\Xt)Jo 



By applying Lemma 16.81 with p{x) = sinh (x), we obtain 

Xt rXt ^ rXt 

{ax + b)'^ smh^ (x) dx < 2 / u(a;)^ sinh^(x) + K{t) / u{xY dx 

Jo Jo 

where K{t) ~ K{t) ■ sinh^(A'i). By combining this with (HOI) and ([^T|) we have 

Xt or^ -V- pXt pXt 

u^{x) dx < - — / u^{x)smh^{x) dx + W-XfK{t) / u^{x) dx. 



sinh"(Xt) Jo 

Arguing as above, one sees that Xt ■ K{t) tends to zero as t tends to zero. It follows 
that for all sufficiently small i, we have 

dx < ^ — — f (x) sinh^ (x) dx. 

- sinh2(X0 Jo ' ' ^ ' 

Note that CXt/ sinh^(X4) < t" for sufficiently smaU t. 
Since 

inf sinh2(x) = smh^iXt) = t.-"^' - 1, 

[Xt.Yt] 

we also have 

rYt ^2e pYt 

u^(x) dx < — / u^(x) sinh^(x)dx, 

'Xt l-t " Jxt 

The claim follows. □ 



7. Generic simplicity for polygons 

In this section we combine the convergence of analytic eigenvalue branches with 
the convergence of eigenfunctions to generalize our earlier results |HlrJdg07| on 
spectral simplicity of simply-connected polygons to several settings. 
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7.1. Slits and simplicity. Our results on spectral simplicity for polygons depend 
on the following. 

Proposition 7.1. Let T.t be a slit of length 2t centered at a point p belonging to 
the interior of a Lipschitz domain C M^. Let D be a measurable subset of dfl 
and let D[ C fls be either D or D UY.. If the spectrum of q on Hj~){fl) is simple, 
then for all but countably many t, the spectrum of q on H^i{fl^^) is simple. 

Proof. By Thcorcni l4.21 the eigenvalues vary analytically for t > 0. Hence it suffices 
to show that there docs not exist a real-analytic eigenvalue branch Et such that 
the dimension of the associated eigenspace Vt is greater than 1 for each t > 0. 

Suppose that ipt and ip^ € Vt denote normalized real-analytic eigenfunction 
branches that are mutually orthogonal for each t. By Theorem 16.51 the corre- 
sponding eigenvalue branch converges to some Eq when t goes to 0. By Theorem 
13. 1[ ipt and converge to eigenfunctions ip and '0* on fls with the same eigenvalue 
Eq. Since ipt and V't are normalized and orthogonal for each t, the limits ip and 
ip* are orthogonal. But this contradicts the assumption that the spectrum of fls is 
simple. □ 

Remark 7.2. The proof applies equally well if 51 is a slit domain with simple spec- 
trum. Thus one can iterate the procedure, and prove that for all but countably 
choices of slit lengths, the domain slit along finitely many slits has simple spec- 
trum. Actually, if we take all the slits of the same length t, the spectrum is also 
simple for all but countably choices of t but to prove this result one has to adapt 
Theorem 16.51 to the case of several slits. This is easily done using the localization 
principle and the estimates we have proved. 

By combining Theorem 17. 11 Remark 17.21 and the main result of |Hlr JdgOT] , we 
obtain the simplicity of the spectrum of the generic simply connected slit polygon. 
To be precise, let Sn,k be the set of simply connected n-gons P with k disjoint 
slits J^i,Y^2,--' Note that the vertices of P and the cndpoints of each slit 

Si, S2, • • • , Sfc, determine the slit polygon PsiuS2U - uEfc ■ Thus. Sn,k may be nat- 
urally identified with an open subset of K2n+2fe+2fe^ j^-^ particular, Sn,k inherits a 
natural affine structure and Borel measure. 

Theorem 7.3. If n > A and k > 0, then almost every slit polygon in Sn,k has 
simple Dirichlet or Neumann spectrum. 

Proof. Let P be a simply connected n-gon. Let C = {ci,...Ck} be a set of k 
distinct points in P, and let L = {£1, . . . ,ik} be a set of (not necessarily distinct) 
lines that pass through the origin in IR-^. Let A{P, C, L) be the set of slit polygons 
-PEiuE2U- -uEfc where is centered at Ci and is parallel to £i. Note that a point in 
y4.(P, C, L) is determined by the lengths of the slits. 

The sets A(P, C, L) provide a natural smooth foliation of Sn.k by fc-dimensional 
planar sets. By Proposition 17.11 and Remark 17.21 if P has simple spectrum, then 
almost every slit polygon in A{P, C, L) has simple spectrum. By the main result of 
IHlrJdgOf] , almost every simply connected n-gon has simple spectrum. The result 
follows by integrating transversely to this foliation using Fubini's theorem. □ 

7.2. Polygons with mixed boundary conditions. Let P be a simply connected 
polygon, and let D be a union of a set of edges in the boundary, dP, of P. Recall 
from ^that the eigenfunctions of the Dirichlet energy q on II^{P) satisfy Dirichlet 
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conditions on D and Neumann conditions on dP\D. Wc say that (P, D) has simple 
spectrum if and only if the quadratic form q on {P) has simple spectrum. 

Let vi,V2-, ■ ■ ■ ,Vn be a cyclic ordering of the edges of a simply connected polygon 
P. Let Ci denote the boundary edge joining the vertex Vi to Wi+iQ Given a subset 
(T c {l,2,...,n|, let 'Pn,a denote the set of pairs {P, D) where P is a simply 
connected n-gon|f| and 

D = \Je,. 

The set Pn,cr can be naturally identified with an open subset of K^" and hence has 
a natural Borel measure and afHne structure. 

Theorem 7.4. // n > 4, then almost every polygon in Vn,a has simple spectrum. 
Proof. By arguing as in HlrJdgOT] , we see that it suffices to construct one polygon 



in Vn.a that has simple spectrum. 

We first construct such a polygon in the case of alternating boundary conditions, 
that is, we suppose that i G a if and only if i + 1 ^ ct where as usual n + 1 = 1. 
Note that n is even in this case. 

Let Q be the rectangle [0, a] x [0,1], and let D = {0,1} x [0, a]. That is, we 
consider the eigenvalue problem on Q with Dirichlet conditions on the vertical 
edges and Neumann conditions on the horizontal edges. The set of functions of the 
form 



sm 



/miTT \ 

' ' ("^27r • y) 



where mi > and 7712 > are integers, is an orthonornal basis of eigcnfunctions. 
Thus, the set of eigenvalues is {(mi7r/a)^ + (TO27r)^}. In particular, {Q,D) has 
simple spectrum if and only if ^ Q. 

Let ci, C2, . . . Cfc be fc = ?i/2 — 2 distinct points on the horizontal segment [0, a] x 
{1/2}. By Theorem 1 7. 1[ there exist horizontal slits Si, . . . , centered respectively 
at Ci such that (Qsiu - uSi- , D') has simple spectrum where D' = 13 U Ei U • • • U S^. 
See Figure [H 

The slit domain Qsiu - uSfc has a reflection symmetry t induced by {x,y) 1— s- 
(a;, 1—y). A standard argument shows that restricting functions on Qsiu-'-u^k to the 
rectangle Q' ~ [0, a] x [0, 1/2] defines a bijection between r*-invariant eigenfunctions 
of q on ^^^(QsiU - uEfc, i^') and the eigenfunctions of q on H]j,{Q'). 

Therefore, since the spectrum of (QsiU - uEfc , is simple, the spectrum of 
(Q', D') is also simple. By placing vertices at the endpoints of the slit, we may 
regard the rectangle Q' as a polygon with n = 2fc + 4 vertices. From this viewpoint, 
the boundary conditions given by D' alternate as desired. 

Given a general subset a C {1,2,...,??}, let F be the finite set obtained by 
identifying i and i + 1 if they either both belong to a or if they both do not belong 
to (T. The cyclic ordering of {1, 2, . . . , r?} induces a cyclic ordering on P, and hence 
we may identify F with {1,2, . . .r?'} for some n' < ??. Note that n' is either 1 or 
even. Let a' C {1, 2, ... r?'} denote the set obtained by identifying i G a with i + 1 
if it also belongs to cr. 

The set a' C {1,2,..., ??'} is alternating, and hence, if n' > 4, there exists an 
element {P' , D) in T'n',a' with simple spectrum. By judiciously adding vertices to 



^Here Sn joins Vn and vi 
pi 

many line segments and whose interior is simply connected 



'^By simply connected polygon, we mean a compact set whose boundary consists of finitely 
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(a,l) 
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Si 



S2 



(0,0) 



(a,l) 



Figure 1 . The slit rectangle 



the boundary edges of the n'-gon P, we obtain an n-gon P' such that the set D 
corresponds to a. Thus {P, D) is an element of Vn,a with simple spectrum. 

If n' = 2, then we consider the eigenvalue problem on {Q, D) where Q = [0, a] x 
[0, 1] and D is either {0} x [0, 1] or its complement. In either case a basis of 
eigenfunctions can be constructed by taking products of sines and cosines. By 
making an appropriate choice of a, one finds that the spectrum of (Q, D) is simple. 
By adding vertices appropriately to the boundary edges of Q, we obtain (Q', D) e 
'Pn,a with simple spectrum. 

The case n' = 1 is the main result of [HlrJdgOT] and follows from the same 
construction. □ 

7.3. Mutliply connected polygons. We first make precise the definition of mul- 
tiply connected polygon. Let Pq, Pi, . . . , Pk be a finite collection of simply con- 
nected polygons. Assume that 

(1) Pi is contained in the interior of Pq for all i > 0, and 

(2) P^ nPj =0 for aU i,j >0, j, 

The multiply connected polygon P determined by Pq., Pi . . . , Pk is obtained by re- 
moving the interiors of the polygons Pi, ... from the polygon Pq. In other words, 
P is the closure of 



Let ft = (no,ni,-- - ,nk) denote a vector of integers with Ui > 3 for each i. 
Let 'P(n) denote the set of all collections of polygons Po, Pi, . . . , Pfc satisfying (1) 
and (2) above and such that for each i, the polygon Pi has Ui (ordered) vertices 
Vi^i, . . . jVi^m- Since the (ordered) vertices determine the polygon. P(n) is naturally 
in bijective correspondence with an open subset of where d — 2nQ + • • • + 2nk- 
In particular, V{n) inherits an affine structure and a Borel measure. 




Proposition 7.5. The space V{fi) is path connected. 
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Proof. In |Hlr JdgOT] , we proved that the space Vn of simply connected n-gons is 
connected using a construction which we will call 'deleting a vertex'. In particular, 
given a simply connected n-gon P, we constructed a linear path 1 > P{t) of n-gons 
in Pn with P(0) = P and such that P(l) has three consecutive vertices that belong 
to the same boundary edge. See Figure [2l We can regard the polygon P(l) as 
an element of Vn-i- Thus, since the space of triangles, V^, is connected, Vn is 
connected for n > 3 by induction. 



V2 




Figure 2 . Deleting a vertex 

Suppose now that P e 'P{n) is defined by the simply connected polygons Pq, Pi, 
. . . ,Pk where > 3 for some i > Q. Let 1 1-^ Pi{t) denote the path from Pi — Pi(0) 
to a polygon Pi(l) having three consecutive vertices on the same edge. Let 1 P{t) 
denote the path of polygons determined by Pi, ... , Pi_i, Pi(i), Pi+i, • • • , Pfc- The 
polygon P(l) may be regarded as an element of Vino, . . . , rti — 1, . . . , n^). 

Therefore, by induction on the vector ft, we find that it suffices to prove that 
the space 'P{n,3,3,3, . . . ,3) is connected. By rescaling the 'interior' triangles if 
necessary, we may delete vertices of the polygon Pp to obtain a path to an element 
of 7^(3, 3, 3, 3, ... , 3). An elementary argument then gives that 7-'(3, 3, 3, 3, . . . , 3) is 
connected. □ 

Theorem 7.6. If uq > 4, then almost every polygon in T'{n) has simple spectrum. 

Proof. Since V{n) is connected, by arguing as in |Hlr JdgOf] , we see that it suffices 
to construct one polygon in V{n) that has simple spectrum. 

Let P belong to 7^(4, 3, 3, 3, . . . , 3). By judiciously adding vertices to the bound- 
ary edges of P, we may regard P as an element of V{no, ni, n2, n^, . . . , Uk) where 
no > 4 and > 3 for i > 0. Thus, it will suffice to prove that there exists some 
P £ P(4, 3, 3, 3, . . . , 3) such that P has simple spectrum. 

By Theorem 17.31 there exists a quadrilateral Q with k slits Ei, S2, • ■ • so 
that (5EiuE2U - uSfc has simple spectrum. For each i, choose a point pi so that 
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the convex hulls, Pi, of {pijljSi satisfy conditions (1) and (2) above. Let rrii 
be the midpoint of E.; and define the path Xi{t) = tpi + (1 — t)'mi. Define Pi{t) 
to be the convex hull of {xijljEi. Let P{t) be the multiply connected polygon 
defined by Q, Pi{t), P2{t), . . . , Pkit). See Figure [3 Note that P{t) is an element of 
7^(4, 3, 3, . . .,3) for t > and that P(0) corresponds to ((5,I]i,I]2, . . . ,S„). 




Figure 3. Perturbing a slit into a triangle 

We define a linear path t t—^ ft oi piecewise linear homeomorphisms mapping 
P(l/2) onto P{t) for each t S [0, 1/2]. In particular, let Qi{t) be the quadrilateral 
Pi{l)\Pi{t), and let /( be the 'obvious' piecewise linear homeomorphism that maps 
g,(l/2) onto Qi{t). Define /* : P(l/2) ^ P{t) by 

..^N / ^, ifx^UQ.(l/2) 
' ~ I fti^), ifxGg.(l/2). 

By Lemma 2.1 in |Hlr Jdg07| , the eigenvalues of P{t) vary real-analytically in t. 
Since the spectrum of P(0) is simple, there exists to > such that the spectrum of 
P{to) is simple. The claim is proven. □ 
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Appendix A. Mathieu functions 

For the convenience of the reader we provide basic facts about Mathieu functions 
that are used in the present paper. For additional information, we refer the reader 
to |Mth68| and [MrsFsh] . Our approach is based on analytic perturbation theory. 
See Example VII.3.4 in [Katoj . 

Suppose that A?/; = E ■ ip. A straightforward computation shows that 

(42) -(92 +d^) [i^oFt) = f ■E-{cosh{zf - cos^ (6*) o F*). 

Let h = t\fE. The method of separation of variables leads one to consider the 
operator 

Ah = + h^cos'ie) 

acting self-adjointly on i^(R/(27rZ), 

Proposition A.l (Angular Mathieu functions). For each z € Z, there exist unique 
real- analytic paths Vi : R — > L^(M/ (27rZ), d6) and &i : R — > M so that for each ft, g R 

(a) Vi{h) is an eigenfunction for Ah with eigenvalue bi{h): 

Ah{vi{h)) = b^{h) -Viih). 

(b) Vi{h) has unit norm in i^(R/(27rZ), rffl), 

(c) the span of {vi{h) \ i = 0,1,2...} is dense in L^{R/{2TTZ),dd), 

(d) For i odd 

v^{h){e) ^ 7r-5 . sin(i6i) + 0(/i^), 

and for i even 

v,{h){e) = 7r"5 . cos(i6l) + 0(h^). 

(e) For i odd 

b,{h) ^ i^ + 2/i2 + 0{h'^), 

and for i even 

h{h) ^ + ^h^ + 0(ft4) 

(f) For any i and any h, we have hi[h) > 6^(0). 

Remark A. 2. The functions Vi{h) arc known classically as angular Mathieu func- 
tions, and we will use indifferently the notations Vi{h) and Vi_h- See for example 
jMrsFshj . 

Proof. The path h s- is an analytic family of compactly resolved operators. 
Hence the existence of paths Vi and hi satisfying (a), (b), and (c) follows from 
Theorem VII.3.9 on page 392 in [Kato] . 

The fc^-eigenspace of Aq = dg is spanned by Bi = {sm{i9) , cos{i9)} . Let Pi be 
the orthogonal projection onto this eigenspace. Let A denote the second derivative 
of h 1-^ Ah evaluated at ft, = 0. One computes the matrix of PiAPi with respect to 
the orthonormal basis Bi to be 

(f l)^ 
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Uniqueness of the paths as well as property (d) then follow from analytic pertur- 
bation theory [Katoj . The derivative of bi is given by 

Thus hi is increasing for positive h and decreasing for negative h showing that h = 
is a global minimum. □ 

Corollary A. 3. Let ip be an eigenvector of A with Dirichlet boundary condition 
on the slit (resp. Neumann) then 

i 

where 

Viji = Vi{h) 

are as above, and u = Ui_h satisfies the ordinary differential equation 
(43) - u"{z) + ( - cofil? [z)) -uiz) = 0, 

with Dirichlet boundary condition at z ~ (resp. Neumann). 
Proof. Each eigenfunction -0 is smooth and hence by Proposition lA.il 

tl;{z,0) = y^M,jt(z) ■ Vi^hjd). 

i 

where h'^ = t'^E, and 

u,,h{r) - / ^{z,0)-v,,h{O) dO. 
Jo 

One obtains ^5)) by using and integrating by parts. □ 
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